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MINIMUM-MASS ISOTROPIC SHELLS OF REVOLUTION SUBJECTED 
TO UNIFORM PRESSURE AND AXIAL LOAD* 

By W. Jefferson Stroud 
Langley Research Center 

SUMMARY 

A method is presented for calculating the shape and thickness distribution of a 
minimum-mass shell of revolution for given coordinate and slope end conditions and given 
values of axial load and pressure load. The method makes use of linear membrane theory 
to relate the loading and shape to stress resultants and uses the strain energy of distortion 
(Von Mises) yield condition to relate the stress resultants to the thickness. The mass can 
then be expressed as an integral which is a function of shape. The integral is minimized 
by a Ritz procedure together with mathematical programing methods — that is, the shell 
shape is written in the form of a series of functions with undetermined coefficients, and 
the mass is minimized with respect to these coefficients. The shape, thickness distribu- 
tion, and stress distributions are given for several examples. 

INTRODUCTION 

In the design of aerospace structures, it is important to keep the structural mass 
low while providing adequate structural integrity. Efforts to develop low-mass, high- 
strength aerospace structures have produced many efficient structural concepts, including 
stiffened- skin construction, sandwich construction, and composite materials. For certain 
structural applications the unstiffened isotropic shell of revolution is efficient. The pres- 
ent report examines unstiffened isotropic shells of revolution designed for minimum mass 
by adjusting both the shell shape and the thickness distribution. The loading is assumed 
to be a combination of axial load and uniform hydrostatic pressure. 

Studies of low-mass isotropic shells of revolution have been performed from various 
viewpoints. The shell middle surface has been assumed to be defined and optimum thick- 
ness distributions of the shell wall for minimum mass have been determined in refer- 
ences 1 to 9. The shell wall thickness has been assumed constant and the optimum 
middle -surface shape has been determined in references 10 and 11. Finally, minimum- 
mass shells for which the wall thickness and middle -surface shapes are optimized 

♦Part of the information presented herein was included in a thesis entitled "The Use 
of a Ritz Procedure in Determining Minimum Weight Shells of Revolution," submitted in 
partial fulfillment of the requirements for the degree of Doctor of Philosophy in 
Engineering Mechanics, Virginia Pol 3 h:echnic Institute, Blacksburg, Virginia, June 1967. 



concurrently are considered in references 12 to 14. In all this previous work the failure 
criterion is assumed to be material yielding; buckling or wrinkling of the shell wall is not 
taken into consideration. Furthermore, except in reference 11, the shell middle -surface 
shapes are limited to very specific classes such as ellipsoidal or cylindrical, or to shapes 
characterized by only one or two shape parameters. Reference 11 considers shells which 
are closed at one end but which have a general middle -surface shape. 

The object of the present report is to give a technique for calculating both the shape 
and the wall thickness distribution of a minimum-mass isotropic shell of revolution for 
given allowable geometric end conditions and given values of axial load and pressure load. 
This work is related, therefore, to that in references 12 to 14; however, the present work 
allows very general middle-surface shapes rather than being restricted to specific classes 
of shapes. In the present report, as in the previous work, the failure criterion is assumed 
to be yielding of the material, and buckling of the shell wall is not considered. Strictly 
speaking, since buckling is not considered, the present work is limited to thin-walled 
shells which are in tension at every point in the middle surface. In addition, variations 
in enclosed volume are not considered in the derivation of the minimum-mass configura- 
tions. In many cases, inequality constraints can be imposed during the synthesis process 
so that the resulting design will have no compressive stresses and/or so that the resulting 
design will enclose at least a minimum volume. Designs with inequality constraints on 
the stresses are considered in this report, but constraints on volume are not considered. 

Linear membrane theory is used to derive all shapes and thickness distributions. 
Results for several examples are presented by means of nondimensional curves and 
tables. 


SYMBOLS 

The physical quantities in this paper are given both in the U.S. Customary Units 
and in the International System of Units (SI) , Factors relating the two systems are given 
in reference 15; those factors used in the present paper are given in the appendix. 

aj^,a 2 ,a 3 constants in expression for shell shape which are used to satisfy the 
boundary conditions (see eq. (32)) 

b depth of pressure-vessel head 

b dimensionless depth of pressure-vessel head, b/y^ 
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constants in expression for shell shape which are used to minimize 
the mass (see eqs. (31), (32), and (33)) 



m 


m 


Ne, 

N0 

P 

Rl, R2 


Xi 


axial load per unit circumferential length applied to the shell at x = 0, 
positive when tensile 

mass of shell 

dimensionless mass, mo^^ppy^^ 

stress resultants in circumferential and meridional directions, 
respectively 

dimensionless circumferential stress resultant, 2N^^pyQ 
dimensionless meridional stress resultant, 2N^y^py^ 
pressure, positive when internal 
shell radii of curvature (see eqs. (10) and (11)) 
distance along meridian 

dimensionless distance along meridian, s jy^ 
thickness of shell wall 
dimensionless thickness, 2tCT(./pyQ 

coordinate along shell axis of rotation, with origin at left end of shell 
value of X at right end of shell 
dimensionless axial coordinate, x/y^ 
dimensionless length of shell, Xj^y^ 

shell radius; radial distance from shell axis of rotation to point on shell 
surface 


t dy 
^ "a 

dx 
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value of y at x = 0 (see fig. 4) 



yo 


y 



y'* 



y 


o 


p 


%,o 

^ 4 > 

<t> 


CO 


dimensionless shell radius, Y/Yq 


slope of meridian at x = 0 
mass density of shell wall material 
critical stress 
meridional stress at x = 0 
meridional stress 
circumferential stress 

angle between tangent to shell meridian and a line parallel to shell axis 
(see fig. 4) 

value of <p at x = 0 
loading parameter, 2Fo/pyQ 

DESCRIPTION OF PROBLEM 


The minimum-mass design technique discussed in this report is applied to two 
classes of shell shapes: the transition section and the pressure-vessel head. The first 
shape to be examined is the transition section. 
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As used herein, a transition section is a shell of revolution which is used to join 
two other shells of revolution in such a way that all three shells have the same axis of 
revolution. An example is shown in figure 1. If the axis of revolution is identified as 
the X-axis, then the transition section extends from x = 0 to x = Xj. The transition 
section can be designed to support axisymmetric loadings such as normal pressure and 
axial loads. With respect to this first class of shell shapes, the problem is to determine 
the shape and thickness distribution for a minimum-mass transition section designed to 
support a given combination of pressure and axial load. 



Figure 1.- Transition section joining a cylinder on the 
left vith a sphere or cone on the right. 

The second shape to be examined is the pressure -vessel head. The head is used 
as an end closure for an arbitrary shell of revolution. An example is shown in figure 2. 
Here the problem is to determine the shape and thickness distribution for a minimum- 
mass pressure -vessel head. This problem is complementary to the minimum-mass 
transition problem because combinations of minimum-mass transition sections and 
minimum-mass pressure-vessel heads provide complete, closed, minimum-mass shells. 



' — Conical pressure vessel 

Figure 2.- Pres sure -vessel head used as closure for 
a conical or ellipsiodal pressure vessel. 
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ANALYSIS 


The mass m of a thin shell of revolution (see fig. 3) is given by 

m = 2pw r ^ ty)ll + (y’)^ dx (1) 

'^x=0 

in which p is the density of the material (assumed constant), t is the thickness, and 
y(x) is the curve which is rotated about the X-axis to generate the surface. It is 



Figure 5.- General shell of revolution hounded hy 
planes x = 0 and x = x^. 

necessary to express the thickness as a function of y(x) and its derivatives and the type 
and intensity of loading. Such a relationship can be established by means of a critical 
stress condition and shell theory. 


Critical Stress Condition 

In this report it is assumed that the strain energy of distortion (Von Mises) yield 
condition is an adequate failure criterion. Without intending to dismiss as unimportant 
the problem of instability, buckling is not considered herein. The stresses at each point 
of the shell are thus related by 

ai2 - aj02 + (2) 

in which and 02 are principal stresses and is some experimentally deter- 
mined critical stress such as the yield stress. From linear membrane shell theory, if 
an axisymmetric stress condition is assumed, the principal stresses are given by 


a 


1 " 


t 



(3) 
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so that 


t = 



N()N0 


+ n/) 


1/2 


( 4 ) 


Loading and Stress Resultants 


Consider an axisymmetric shell acted upon by a uniform hydrostatic pressure 
loading p and by some externally applied axial load per unit circumferential length, 
Fo, applied at x = 0 as shown in figure 4. Equilibrium of forces in the axial direction 
requires that the meridional stress be given by 


CT^t27ry cos (p = cos (p^ + 2irp y dy 


(5) 


in which the subscript o is associated with points on the circumference of the shell at 
X = 0, and all other quantities are shown in figure 4. 



Figure 4.- Axisyinmetric loading acting on shell of revolution. 

The term ^^hydrostatic pressure” means that the pressure is applied to all surfaces 
of the shell, including the ends. It is assumed that pressure-vessel heads at axial loca- 
tions beyond x = 0 and x = x^ transmit the axial pressure load into the shell. The 
axial load per unit length of circumference, Fq, is an additional, externally applied load 
and is not directly associated with the pressure. 

With this definition of hydrostatic pressure, it follows that the meridional stress 
o appears in equation (5) is given by 
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^(p,o^o 


F 

° 2 /cos 00 


A parameter w which describes the loading is defined as 


CO 


= '^o ■ ^ ' py^ 


( 6 ) 


(7) 


The loading parameter co provides a convenient method for nondimensionalizing with 
respect to the pressure. It is the ratio of the externally applied axial load at x = 0 to 
the axial component of the pressure load at that same cross section. Since the pres- 
sure p can be internal (positive) or external (negative) and since Fq can be tensile 
(positive), zero, or compressive (negative), the parameter co can be positive, zero, or 
negative. In defining the loading parameter co in such a way that the pressure p is 
in the denominator, it is assumed that p is not zero. The case of p = 0 is handled 
separately as a special case. 


From equations (5), (6), and (7), the stress can be expressed as 


V = = 


PYr 


2y 


\CjO + 



( 8 ) 


The equation of equilibrium of forces normal to the surface of the shell is 


Rl ^R2 


P 


(9) 


By making use of equations (8) and (9) and noting that 

,3/2 


Rl = - 




and 


R2 = y]li + (y')^ 


an expression for Ng can be formed: 

2 / 


Ne = py\/l + (y’)^ + 


pyr 


CO + 


yoy\li+(ff 

The following nondimensional barred quantities are introduced: 

y = -^ 
yo 

x = iL 


(10) 


( 11 ) 


( 12 ) 


(13) 

(14) 
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Then 


y' = y 

(15) 

v" = y_ 
^ yn 

(16) 


in which the primes denote differentiation with respect to x and the dots denote differ- 
entiation with respect to x. 

Nondimensional stress resultants and Ng can be expressed in terms of the 

loading parameter co and the nondimensional shape as follows: 






<P py. 


(17) 


A nondimensional thickness t can be expressed as 

2ta /— o — 2\1/2 

t = -NgN^+N^ ) 


(18) 


(19) 


By using equation (1), with equation (19) relating the thickness to the stress resul- 
tants and equations (17) and (18) relating the stress resultants to the geometry and 
loading, the nondimensional mass m can be expressed as 


ma,^ 


m s 


r^=^i -_,| 71 ^ _ 

" J- n 


8 


( 20 ) 


TrpPYo" 

In the special case for which the pressure p is zero, the following equations apply: 


N0 = 




(21) 


= ( 22 ) 

l|i +(y) 

The thickness is calculated from equation (4) and the mass is calculated from equation (1). 

As formulated in equation (20) the shell wall thickness varies along the meridian. 
Minimum-mass shells of uniform thickness can be generated by basing the thickness upon 
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the average stress condition in the shell or upon the most critical stress condition in the 
shell. (See, for example, refs. 10 and 11.) The present report considers only those 
shells for which the thickness varies along the meridian. 

A method for obtaining the shell shape y(x) which minimizes m is presented 
after the required boundary conditions are examined. 


Boundary Conditions 

The number and types of boundary conditions which are appropriate to this opti- 
mization problem can be obtained by use of the calculus of variations. Let equation (20) 
be written as 


_ (- -• - 

m \ fiy,y ,y jdx 


x=0 

where, from equations (17) to (20), f is given by 


(23) 


_ 1 (c.-ny^f[i4-(yf] _ _ 2 L + y2) H + (ff\ ^ + y^)^(y ') 


a /2 


4yy (o) + y^) + 4y^ J + (y ) ^ / y]ji + (y ) 


(24) 


The variational equation for the functional m 

pX=Xl 


is of the form 


6m = 


|)ay . 


.9y> 


8y 


dx = 0 


(25) 


In order to examine the types of boundary conditions that are required, the second term 
in the integrand is integrated by parts once and the third term is integrated by parts 
twice to give 


6m = 




+ 



- X=Xj 

Sy* 

-^x=o 



9y 




6y dx = 0 


(26) 
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The boundary conditions which are given in equation (26) are: 
At X = 0, 


af 





af 


9y 


= 0 


At X = Xj, 

af 





af 




= 0 


or y prescribed 
or y prescribed 


or y prescribed 
or y prescribed 


(27) 

(28) 


(29) 

(30) 


These are all the legitimate choices for boundary conditions. The boundary conditions 
that are used in the examples presented in this report are the prescription of y and y 
at each end. 

For the case of o) = -1.0, examination of equation (24) shows that di/dy is zero 

at x=0 ^since yQ = ^=l.(^. Therefore, from equation (28), the slope y cannot be 

prescribed at that end. Note that when o> = -1.0 the externally applied axial load at 
x = 0 cancels the axial component of the pressure load at the same end so that is 

zero there. Thus, equilibrium can be maintained even if a slope discontinuity exists. 

It is emphasized here that deflections and deformations are not considered in this 
report. The coordinate and slope boundary conditions refer to the shell shape at the ends 
and do not refer to deflections away from some initial shape. 


MINIMIZATION OF MASS BY THE RITZ PROCEDURE AND 
NONLINEAR PROGRAMING METHODS 


The problem of finding the shape and thickness distribution of a minimum-mass 
shell of revolution has been cast in the form of finding the function y(x) which mini- 
mizes a functional of the form ^ f(y,y*,y")dx. In this report Ritz's method, together 
with nonlinear programing techniques, is used in the minimization. 

Two types of shell shapes are considered. The first t 5 q)e is denoted transition sec- 
tion. For this type of shape, the shell radius y is not zero at either end. The second 
t 5 q>e of shell shape is the pressure -vessel head. 
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Transition Sections 


The examples for transition sections presented in this report include only shapes 
that are symmetric about the midlength Xj/2. An equation which can be used to describe 
any smooth transition section of this type is 



n=l,2,3 


c„ sin2 ^ 
n xj 


(31) 


The first two terms on the right-hand side of equation (31) are used to satisfy the required 
coordinate and slope boundary conditions. The factor y^ is the slope at x = 0. The 
constants c^ are used to minimize the mass. The trigonometric series can describe 
any smooth curve that is symmetric about the midlength and that has zero coordinate and 
slope boundary conditions. Regardless of the values of Cj^ the trigonometric series 
makes no contribution to the coordinate and slope boundary conditions that are applied to 
the transition section. 


Although the transition sections presented in this report are limited to shapes that 
are symmetric about the midlength, a more general expression which could be used to 
generate unsymmetric shapes is 


y = 1 + a.x + a„x 


-2 


+ a3X 


n=0,l,2 L_ 


cos 


HIE _ cog (n j-2)7ri^ 


^1 


(32) 


The constants a^, 3 . 2 , and a^ are used to satisfy the boundary conditions. 


Pressure-Vessel Heads 


The pressure-vessel heads considered herein act as closures for axisymmetric but 
otherwise arbitrarily shaped pressure vessels. At the junction of the head and vessel 
the coordinates and slope of the pressure-vessel head are assumed to match those of the 
pressure vessel. All loads acting on the pressure-vessel head are assumed to result 
from uniform hydrostatic pressure; therefore the loading parameter co is zero. But 
this requirement on 00 does not extend to the pressure vessel to which the head is 
attached. The pressure vessel can carry an externally applied axial load, as shown in 
figure 5. 


Pressure vessel, 




j Externally applied 

axial load 

Pressure-vessel head, o = 0 
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Figure 5*-* Pressure-vessel head attached to cylindrical pressure 
vessel ■which is carrying an externally applied axial load. 



An expression for y(x) which can satisfy the boundary conditions and which forms 
a smooth closure for a pressure vessel may be written as 



in which b is the dimensionless depth of the head. Two of the arbitrary coefficients 
are used to satisfy the coordinate and slope boundary conditions at x = 0, the junction of 
the head and vessel. The remaining coefficients and the depth of the head are used to 
minimize the mass. 


The series given in equation (33) was chosen because it provides smooth stresses 
at X = b, the apex of a pressure-vessel head. Mathematically the smooth-stress 
criterion is given by 


dN. 

lim 

- - ds 
X— b 


= 0 


(34) 


and 


^d> 

lim = 0 

x^h ds 

where s is the dimensionless distance along the meridian. 

Equations (34) and (35) guarantee that the thickness is also smooth at the apex 
that is, 

lim ^ = 0 
x^h ds 


(35) 


(36) 


Nonlinear Programing Methods 

The values of Cj^ in equations (31) and (33) were obtained by nonlinear programing 
methods. The c^ are the design variables and the mass in is the objective function. 
The direction-generating algorithms that were used include random, gradient, sectioning, 
Newton-Raphson, and Fletcher-Powell. After the direction-generating algorithm was 
used to calculate a set of direction cosines for a move, a one -dimensional search was 
carried out to find the minimum-mass design in that direction. 


In the random method the direction cosines are generated from random numbers. 

In the gradient method the move direction is the negative gradient direction. In the sec- 
tioning method the design variables are improved sequentially rather than simultaneously. 
As used in this report, the Newton-Rapjison method (ref. 16) searches for the Cj^ that 
are solutions to the set of equations = 0. The Fletcher-Powell method is described 


in references 17 to 20. 


acr 
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The computer program was written in such a way that several direction-generating 
methods could be used in a single computer run. Several methods were used to start 
designs. The Fletcher-Powell method was used to complete all designs having 10 or 
more design variables Cjj. 


RESULTS AND DISCUSSION 
Transition Sections 

Linear membrane shell theory, the Von Mises yield condition, and equation (31) 
were used to generate typical minimum-mass shapes and thickness distributions for 
transition sections having four combinations of geometric end conditions, all of which 
are symmetric with respect to the midlength. One transition section has a length equal 
to the radius at each end and has zero-slope end conditions. Another transition section 
has a length equal to 1.5 times the radius at each end and has zero-slope end conditions. 
A third transition section has a length equal to 1.5 times the radius at each end; the slope 
at the left end is 0.5 and the slope at the right end is -0.5. The fourth transition section 
is a special case (w = -1.0) for which the slopes at the ends cannot be prescribed in 
advance. For this fourth transition section, the length is 1.5 times the radius at each 
end. Results are shown in figures 6 to 19. 

The discussion of transition sections includes, first, a presentation of shapes, 
thickness distributions, and stress resultants. Next, certain aspects of these results 
are examined. Finally, an indication of potential mass savings is presented. 

Shapes, thickness distributions, and stress resultants.- Three slope end conditions 
are considered: (1) zero-slope end conditions, (2) nonzero-slope end conditions, and 
(3) the special case co = -1.0 for which the slopes at the ends cannot be prescribed. 

(1) Transition sections with zero-slope end conditions: The transition sections 
shown in figure 6 have a nondimensional length Xj equal to 1.0, which means that the 
length is equal to the radius at each end. The nondimensional meridional shapes for 
several values of o) are shown in figure 6(a). All these shapes neck down in the center 
except for the limiting cases w = ±<» (p = 0), which are cylinders. 

The nondimensional thickness distributions associated with three of the preceding 
shapes, = 0.6, 0.0, and -0.4, are shown in figure 6(b). The thickness for the cases 
CO = ±°o (p = 0) can be calculated from equation (4). The thickness distributions are 
discussed in more detail in a subsequent section entitled "Oscillations in the thickness 
distribution." 

The nondimensional meridional and circumferential stress resultants are shown in 
figures 6(c) and 6(d), respectively. Note that the thickness distributions shown in 
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(a) Nondlme ns zonal meridional shape. 


(c) Nondimens ional meridional stress 
resultant . 



(h) Nondimens ional thickness distribution. (d) Nondimens ional circzimferential 

stress resultant. 

Figure 6.- Nondimens ional meridional shapes^ thickness distributions; and stress resultants 
for several minimum-mass transition sections ^fith = 1.0 and yZ = 0. 





figure 6(b) reflect the stress distributions shown in figure 6(d). The stress resultants 
for u) = ±<x> can be calculated from equations (21) and (22). 

The transition sections shown in figure 7 have a nondimens ional length equal 
to 1.5 and have zero-slope end conditions. This set of transition sections is therefore 
similar to, but 50 percent longer than, the set of figure 6. The increase in length allows 
these minimum-mass designs to neck down in the center more than the previous set. 

The shell radius y at the midlength Xjy/2 of both sets of transition sections is 
shown as a function of the loading parameter u> in figure 8. The upper curve is for the 
transition sections having = 1.0, and the lower curve is for the transition sections 
having Xj = 1.5. For these two sets of transition sections the midlength point is the 
point at which the minimum value of y occurs. Hence, these curves show the minimum 
value of y as a function of cu. The symbols indicate calculated designs. The transi- 
tion sections having Xj = 1.0 are essentially cylindrical for u> > 1.0 and ol> < -3.0. 

The transition sections with x^ =1.5 neck down for a wider range of u> and are essen- 
tially cylindrical for w > 2.0 and u> < -4.0. 

(2) Transition sections with non-zero-slope end conditions: The shapes, thickness 
distributions, and stress distributions for several transition sections having non-zero- 
slope end conditions are shown in figure 9. For these transition sections the slope at the 
left end is 0.5 and the slope at the right end is -0.5. The nondimensional length x^ is 
1.5. Except for the slopes at the ends of the meridional shapes, the results shown in 
figure 9 are similar to the results shown in figures 6 and 7. 



0) 


Figure 8.- Midlength radius 



as a function of loading parameter 


for two types of transition sections with zero- slope end conditions. 


CD 
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(c) Wondimensional meridional stress 
resultant. 


I 1- - . . - --J 

0 .75 1.50 

X 

(a) Wondimensional meridional shape. 




(h) Wondimensional thickness distribution. (d) Wondimensional circumferential stress 

resultant. 


Figure 9.- Wondimensional meridional shapes^ thicimess distributions^ and stress resultants 


for several minimum-mass transition sections -with =1.5 and y 


0 . 5 . 



I 


For transition sections having the same length and end conditions as those of fig- 
ure 9, the shell radius y at the midlength is shown as a function of the loading 

parameter o) in figure 10. The symbols indicate calculated designs. 



for transition sections with X2_ = 1.5 and Jq = 0.5- 


The special case for which the pressure p is zero is presented in figure 11. The 
very great thickness at the ends suggests that a ring may be useful at each end of the 
shell. 

(3) The special case co = -1.0: As was stated previously in the section entitled 
"Boundary Conditions," when the loading parameter cu is equal to -1.0 the slope bound- 
ary condition cannot be prescribed at the left end x = 0. Since the examples for transi- 
tion sections presented in this report include only shapes that are symmetric about the 
midlength, the restriction concerning slope boundary conditions applies to the right end 
X = Xj as well as to the left end. The restriction on prescription of slope boundary con- 
ditions means that in addition to minimizing the mass with respect to the Cn in equa- 
tion (31), the mass must also be minimized with respect to y^. 

An example of a transition section for which cu = -1.0 and x^ = 1.5 is shown in 
figure 12. For this design y^ = -2.784. The design appears to require rings or bands 
at the ends and a relatively uniform wall thickness over the rest of the shell. The mid- 
length radius yfx.^/^ for two cases with w = -1.0 is also included in figures 8 and 10. 


19 







I 
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Figure 12.- Meridional shape ^ thicloiess distribution ^ and stress resultants 
lor minimum-inass transition section with od = -1.0^ = 1*5^ and end 

slopes not prescribed. 

Numerical values of the masses and Ritz coefficients (or design variables) c^^ are 
given, in tables 1(a), 1(b), 1(c), n, and III for the shell shapes presented in figures 6, 7, 9, 
11, and 12, respectively. 

Oscillation^ in the t hickness distribution .- The oscillatory character of some of the 
thickness distributions (for example, see co = -0.4 in fig. 6(b)) arises because a series 
solution is used. For practical construction the oscillatory thickness distributions would 
have to be smoothed in some manner. Part of the problem may arise because the shell 
shape, rather than the shell thickness, is being expressed as a series. The thickness is 
obtained by differentiating this series — a procedure that is known to cause convergence 
problems in some instances. For this reason, it is interesting to see how the thickness 
distribution changes as the number of terms in the shape series is changed. 

Curves which show these changes for the case o) = -0.4 of figure 6(b) are pre- 
sented in figure 13. Because of symmetry, only half the shell is shown. It appears that 
certain trends in the thickness distribution have been established. From x = 0 to 
X = 0.35, the thickness could be smoothed to provide a slowly changing thickness; at 
X ts 0.40 a large peak in the thickness occurs. The mass of the 10-term solution is less 
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than 2 percent more than the mass of the 20-term solution, and the mass of the 15-term 
solution is less than 0.5 percent more than the mass of the 20-term solution. The cor- 
responding changes in the values of y at the midlength of the shell are also small. If 
only the shell shape and the objective function, which is the mass, are considered, the 
series has converged. However, additional terms cause substantial changes in the thick- 
ness distribution. Perhaps including even more terms in the series would cause the 
oscillatory character of the thickness distribution to disappear eventually, but the short- 
wavelength terms which would be added would allow the large peak at x = 0.40 to become 
higher and narrower. In cases such as this, constraints may be imposed to obtain a more 
practical solution. 

Negative (c ompres sive) stresses, buckl ing, and con straints .- In the procedures dis- 
cussed in the preceding pages the critical stress condition used has been the Von Mises 
yield condition, and the possibility of buckling has not been considered even though com- 
pressive stresses may be present in the shell wall. For many loadings and boundary 
conditions it is possible, in designing the shell, to use an inequality constraint which 
eliminates compressive stresses. In other cases, constraints could be used to prevent 
sharp peaks in a stress or thickness distribution and to prevent the shell shape from 
necking down more than a prescribed amount. In this section two examples of con- 
strained solutions are presented. Both are for the co =: -0.4 case of figure 6, and in 
both examples a stress resultant constraint is imposed. 

In the first example, Ng is required to be less than 2.5, The constraint is 
imposed in an attempt to reduce or remove the spikes which occur in the Ng distribu- 
tion at X = 0.4 and x = 0.6. For comparison, Ng = 2.0 for a cylinder. The results 
are presented in figure 14. The constraint turns out to be quite severe. Whereas the 
unconstrained solution is 29 percent lighter than the corresponding cylinder having the 
same length and end radius and subjected to the same loading, the constrained solution is 
only 12 percent lighter than the cylinder. The mass and Ritz coefficients are given in 
table W(a). 

In the second example, H(p is required to be positive (tensile, for internal pres- 
sure) over the total length of the shell. In the unconstrained solution, is negative 

at the center. The results, which were obtained from a 15-term solution, are presented 
as solid lines in figure 15. Portions of the shape, thickness, and stress resultants for 
the 15-term unconstrained design are shown in figure 15 as dashed lines. In this case 
the constraint is not severe. The constrained solution which has no negative stresses 
is 26 percent lighter than the cylinder, whereas the unconstrained solution is 29 percent 
lighter than the cylinder. The mass and Ritz coefficients are given in table IV(b). 
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(a) Nondimensional meridional shape (h) Nondimen si onal stress resultants, 

and thicloiess distribution. 


Figure l4.- Nondimensional meridional shape ^ thickness distribution^ and stress resultants 
for miniiTium-mass transition section with = 1.0; Yq = ^ ^ -0.4; and the added 

constraint Nq < 2.5* 



(a) Nondimensional meridional shapes (b) Nondimensional stress resultants, 

and thiclmess distribution. 

Figure 15.- Nondimensional meridional shapes; thickness distributions; and stress resul- 
tants for a minimum-mass shell designed to have no negative stresses and for the 
minimum- mass shell having the same loading and boundary conditions but designed with- 
out that constraint. = 1.0; = 0; cd = -0.4. 
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In both of these examples the constraint was imposed by the sequential uncon- 
strained minimization technique (SUMT) method of references 21 and 22, which makes 
use of a penalty function. Additional information on the use of the SUMT method in 
structural optimization problems may be found in references 18, 19, and 20. 

Be nding stresses .- Since bending stresses are neglected in the development of 
these minimum-mass shell shapes and thickness distributions, it is of interest to deter- 
mine whether the stress distributions predicted by bending theory are significantly dif- 
ferent from those predicted by membrane theory. Stresses at the inner and outer sur- 
faces of the shell were therefore calculated by using a computer program based on linear 
bending shell theory (ref. 23) for the case co = 0.2 of figure 7. 

The end radius yo chosen for the shell is 10 in. (25.4 cm) and the length xj is 
therefore 15 in. (38.1 cm). The shell is assumed to be fabricated of a steel which has 
the following properties: modulus of elasticity E = 30 x 10^ ksi (207 GN/m^); yield 
stress Op = 35 ksi (241 MN/m2); Poisson's ratio, 0.3. The boundary conditions used at 
each end of the shell are: no restraint against normal displacement and no restraint 
against meridional rotation. 

The results are presented in figure 16. The stress quantity 
at the inner and outer surface of the shell wall is shown as a function of arc length 
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Figure l6.- The stress 
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at the inner and outer 


surfaces of two minimum-mass transition sections as a function of dis- 
tance s measured along meridian. Yq = 10 in. (25*^ cm); = 15 in. 
(58.1 cm); = 0; to = 0.2. 
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measured along the shell for two values of internal pressure p. The amount by which 
this stress quantity differs from 35 ksi (241 MN/m2) represents the inadequacy of linear 
membrane theory for that loading. The curves show that for low loadings (p = 100 psi 
(0.69 MN/m2)) membrane theory is adequate and that for high loadings (p = 1000 psi 
(6.9 MN/m2)) membrane theory is marginal. The dimensional thickness distribution can 
be calculated by using y^ = 10 in. (25.4 cm) and Uj. = 35 ksi (241 MN/m^) in equa- 
tion (19) together with the thickness distribution t shown in figure 7(b) for co = 0.2. 

For p = 1000 psi (6.9 MN/m2) the thickness at the center is 1.3 in. (3.3 cm) and for 
p = 100 psi (0.69 MN/m2) the thickness at the center is 0.13 in. (0.33 cm). 

The thickness distribution for p = 1000 psi (6.9 MN/m2) was simplified as shown 
in figure 17 and the shell was reanalyzed. The stress quantity at the inner and outer 
surfaces of the modified shell is presented in figure 18 as a function of the arc length 
measured along the shell. The stress quantity remains below the yield stress for the 
entire shell. The dimensionless mass m for the original, unmodified shell is 2.35. 

The thickness simplification increases m to 2.58. 

Mass saving .- An indication of the mass saving that can be realized by using the 
present method is shown in figure 19. For the transition sections shown in figures 6 
and 7 the ratio of the mass m to the mass of uniform-thickness cylinders designed to 
carry the same load is presented in figure 19 as a function of the loading parameter w. 



Figure 17*- Original and simplified thickness distributions, 
p = 1000 psi ( 6.9 MN/m2)j = 1.5; = 0; a> = 0.2. 
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s, cm 



s. in. 

Figure l8,- The stress quantity at the inner and outer surfaces 

of a transition section with simplified thickness distribution as a function of 
distances from left end measured along the meridian, p = 1000 psi (6.9 ilN/m^); 
y^ = 10 in. ( 25 A cm); = I 5 in. (38. 1 cm); y^ =0; cu = 0.2. 



Figure 19- - Ratio of mass m of minimum-mass transition section to mass of 
\ini form- thickness cylinders designed to support same load^ as a function 
of loading parameter cuo Transition sections have zero- slope end 
conditions. 


27 




The radius of the cylinders is y^, the radius of the transition sections at the ends. The 
length of the cylinders is the same as the length of the corresponding transition section. 
The upper curve is for = 1.0 and the lower curve is for = 1.5. The Von Mises 
yield condition is the critical stress condition for the cylinders as well as the transition 
sections. The symbols indicate calculated designs. The maximum mass saving for the 
shorter transition sections, Xj = 1.0, is about 29 percent. The maximum mass saving 
for the longer transition sections, x^^ = 1.5, is about 41 percent. For large positive or 
negative values of w the axial load Fq is the dominant load, and the minimum-mass 
transition section which is symmetric about the midlength and has zero-slope end condi- 
tions approaches a cylinder. 


Pressure-Vessel Heads 

Equation (33) was used to derive the minimum-mass shapes and thickness distribu- 
tions for three pressure-vessel heads. The design parameter that differentiates among 
these three heads is the slope of the head at its junction with the pressure vessel. One 
of the heads is designed for use with a cylindrical pressure vessel, and in tliis case the 




(a) Nondimenslonal meridional shape. (h) Nondimensional thickness distribution. 

Figure 20.- Nondimensional meridional shapes, thickness distributions, and stress resul- 
tants for minimum-mass pressure-vessel heads with slopes at the head-vessel junction 
equal to -0.5> 0, and 0.5. 
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slope of the head at the junction is zero. Each of the other two heads has a nonzero- 
slope end condition and can therefore be used with conical or ellipsoidal pressure vessels. 
One of these heads has a slope equal to 0.5 at the junction, and the other has a slope equal 
to -0.5. 

Shapes .- The shapes for the three pressure-vessel heads are shown L. .igure 20(a). 
As in the case of the transition sections, the X-axis is the axis of revolution. The non- 
dimensional head depths are 0.54, 0.77, and 0.97 for minimum-mass heads having slopes 
at the junction of -0.5, 0, and 0.5, respectively. For these junction slopes, the minimum- 
mass heads obtained by using the series given in equation (33) are 0.01, 0.4, and 5.4 per- 
cent lighter, respectively, than heads obtained by using the best ellipse. 

Thic kness distributions .- The thickness t for each of the three heads is shown as 
a function of the axial coordinate x in figure 20(b) and as a function of the distance s 
measured along the meridian in figure 20(c). Note that dt/ds = 0 at the apex. 



(c) Nondlmensional thiclcness t as a f\mctlon of nondlmenslonal 
distance s along the meridian. 

Figure 20.- Continued. 
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According to the theory considered in this report, a discontinuity in thickness can 
exist between the head and the pressure vessel at their junction. Consider the head which 
has a zero-slope end condition and a cylindrical shell to which the head could be attached. 
If the cylinder is loaded only by pressure and if the thickness of the cylinder, like the 
thickness of the head, is determined by the Von Mises yield condition, then the thickness 
of the head at the junction is about six-tenths the thickness of the cylinder. 

Stress resultants .- Both the circumferential and meridional stress resultants for 
the three heads are shown as functions of the axial coordinate x in figure 20(d). For 
each head the circumferential stress resultant Ng is shown as a solid line and the 
meridional stress resultant is shown as a dashed line. In each case the stress 

resultants are equal, as they should be, at the apex. For the head having an end slope of 
0.5 a large compressive circumferential stress occurs near its junction with the pressure 



(d) Nondlmenslono.1 stress distributions. 
Figure 20, - Concluded. 
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vessel. Compressive stresses are discussed in a previous section entitled ^'Negative 
(compressive) stresses, buckling, and constraints/^ 

The masses and Ritz coefficients for each of these heads are given in table V. 

Discontinuity bending s tresses .- If the slopes of the head and the pressure vessel 
are equal at the junction, then the meridional stress resultant is continuous across the 
junction. But, within the framework of linear membrane shell theory, a discontinuity in 
the circumferential stress resultant will occur if there is a discontinuity in curvature at 
the junction. For the theory considered herein, it is not possible to generate minimum- 
mass pressure-vessel heads so that curvatures match at the junction. The permissible 
boundary conditions are limited to coordinates and slopes. 

In an actual loaded pressure vessel which has a discontinuity in curvature at the 
head-vessel junction, bending occurs at the junction, allowing the stresses and deflections 
to be continuous across the junction. If both the pressure vessel and the head are thin, 
if they are tangent at the junction, and if the discontinuity in curvature is not large, the 
discontinuity bending stresses will be local and the head can be designed on the basis of 
membrane theory together with local reinforcement at the junction. For example, 
according to reference 13, the mass of the reinforcement needed for an ellipsoidal head 
to be used with a cylindrical tank amounts to slightly more than 1 percent of the head 
mass. 

Minimum-mass pressure-vessel heads that are to be used with cylindrical pressure 
vessels are also considered in reference 13. As in the present report, linear membrane 
theory and the Von Mises yield condition are used in reference 13 for design purposes. 
The limitation to coordinate and slope boundary conditions is therefore applicable to the 
work in reference 13. Unlike the present report, which considers unconstrained shapes, 
reference 13 investigates three constrained classes of pressure-vessel heads - ellipsoi- 
dal, torispherical, and Cassinian oval - and minimizes the mass with respect to a single 
shape parameter for each class of shells. In reference 13 the ellipsoidal head is found 
to be the lightest of the three types of heads considered. The present analysis leads to 
the conclusion that for cylindrical pressure vessels (junction slope equal to zero), the 
ellipsoidal head is essentially the lightest head of all possible shapes. 

Cassinian-type pressure-vessel heads are considered in reference 13 because they 
provide zero curvature at the cylinder-head junction and thereby minimize discontinuity 
bending stresses. Two important comments should be made with regard to Cassinian 
heads and minimum-mass pressure-vessel heads. First, within the framework of linear 
membrane theory it is not appropriate to derive minimum-mass pressure-vessel heads 
of the Cassinian type because prescription of curvature boundary conditions is not 
allowed. And, second, since Cassinian heads include only a small portion of the head 
shapes that have a zero-curvature boundary condition, conclusions based on Cassinian 
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heads cannot be extended to all heads that have zero-curvature boundary conditions or to 
other head shapes that are tailored to minimize discontinuity bending stresses. 

To demonstrate both these comments, equation (33) was used to generate a 
pressure-vessel head with the following boundary conditions at the junction of the head 
and vessel: 


y(0) = 1 

(37) 

y (0) = 0 

(38) 

y '(0) = 0 

(39) 


Ten terms were used to define the shape. It is of interest to compare this pressure- 
vessel head with the minimum-mass head that has only equations (37) and (38) as bound- 
ary conditions. The shapes are almost identical. The zero-curvature head weighs about 
3 percent more than the head without this additional constraint. (The lightest Cassinian 
head weighs almost 50 percent more than the minimum-mass head without the additional 
zero-curvature restraint.) The nondimensional circumferential stress resultants for 
both heads are shown in figure 21. The solid curve is for the head that has zero curva- 
ture and slope at the junction. The dashed curve is for the minimum-mass head with 
only the zero-slope boundary condition. Both these curves are based on linear membrane 
theory. For the zero-curvature head, is constrained to be 2.0 at the junction 
X = 0. The stress then falls off rapidly and oscillates about the dashed curve. As more 



Figure 21.- Nondimensional circumferential stress resultants for 
pres sure- vessel head with boundary conditions y(0) =1, 
y*(o)=0, y*'(0)=0 and for minimum- mass head with bound- 

ary conditions y(0) = 1, y* (O) = 0. 
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terms are added to the series, the stress distributions for the two heads will almost 
coincide except near x = 0. In effect, the solution is trying to ignore the improper 
curvature boundary condition. 


CONCLUDING REMARKS 

A method has been presented for calculating the shape and thickness distribution of 
a minimum-mass shell of revolution for given coordinate and slope end conditions and 
given values of axial load and pressure load. Linear membrane shell theory was used to 
relate the loading and shape to the stress resultants, and the strain energy of distortion 
(Von Mises) yield condition was used to relate the stress resultants to the thickness. 
Membrane theory greatly simplifies the calculations and appears to be sufficiently accu- 
rate for the type of shell problems considered. Moreover, the results obtained by using 
membrane theory could serve as the first approximation for calculations based on more 
sophisticated shell theories. 

Two types of shell shapes are considered: transition sections, which are open at 
both ends, and pressure-vessel heads, which are open at one end and closed at the other. 
The results show that minimum-mass transition sections have a tendency to neck down 
away from the ends. They may have large excursions in the thickness distribution. 
Examples which make use of inequality constraints on the stress resultants are included. 
Inequality constraints can also be imposed on the shape and thickness distribution to 
restrain both the necking down of the shape and the large excursions in the thickness 
distributions. When comparisons are made between minimum-mass transition sections 
and constant-thickness cylinders having the same length and radius and designed to carry 
the same load, the minimum-mass transition sections are, for some loadings, 30 to 
40 percent lighter than the corresponding cylinder. The calculations for pressure-vessel 
heads are applicable for use with conical and ellipsoidal as well as cylindrical pressure 
vessels. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 26, 1970. 
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APPENDIX 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


Conversion factors (ref. 15) for the units used in this report are given in the fol- 
lowing table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 
(**) 

Length 

in. 

0.0254 

meters (m) 

Stress 

> ksi 

6.895 X 106 

newtons/meter^ (N/m2) 

Modulus of elasticity 
Pressure 

J 

psi 

6.895 X 103 

newtons/meter2 (N/m2) 


^Multiply value given in U.S. Customary Units by conversion factor to obtain 


equivalent value in SI Units. 

** Prefixes to indicate multiples of units are as follows: 


Prefix 

Multiple 

giga (G) 
mega (M) 

109 

10® 
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TABLE I.- RITZ COEFFICIENTS Cn FOR SEVERAL TYPICAL 
MINIMUM-MASS TRANSITION SECTIONS 


(a) Zero-slope end conditions; Xj = 1.0 




Cn for - 


n 

o 

1 

11 

3 

CO = 0 

CO = 0.6 


m = 1.258 

m = 1.668 

m = 1.823 

1 

-3.961 X 10-1 

-7.570 X 10-2 

-1.095 X 10-2 

2 

1.319 X 10-1 

9.143 X 10-3 

-3.512 X 10-4 

3 

-4.686 X 10-2 

-2.784 X 10-3 

-8.159 X 10-5 

4 

7.312 X 10-3 

9.501 X 10-4 

-2.558 X 10"3 

5 

9.784 X 10-3 

-4.453 X 10-4 

-1.054 X 10-5 

6 

-1.427 X 10-2 

1.973 X 10-4 

-5.948 X 10-5 

7 

1.190 X 10-2 

-1.086 X 10-4 


8 

-6.863 X 10-3 

5.251 X 10-5 


9 

1.891 X 10-3 

-3.079 X 10-5 


10 

1.534 X 10-3 

1.388 X 10-5 


11 

-3.039 X 10-3 

-7.670 X 10-5 


12 

2.960 X 10-3 

1.519 X 10-5 


13 

-1.989 X 10-3 



14 

8.069 X 10-4 



15 

1.132 X 10-4 



16 

-5.837 X 10-4 



17 

6.471 X 10-4 



18 

-4.732 X 10"4 



19 

2.416 X 10-4 



20 

-7.182 X 10-5 
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TABLE I,- RITZ COEFFICIENTS Cn FOR SEVERAL TYPICAL 
MINIMUM-MASS TRANSITION SECTIONS - Continued 


(b) Zero-slope end conditions; = 1.5 



Cji for - 

n 

0) = 0 

u) 0.2 

CO = 0.6 


m = 1.747 

m = 2.353 

ms 2.674 

1 

-6.670 X 10“1 

-2.345 X 10-4 

-6.012 X 10-2 

2 

2.587 X 10-1 

4.199 X 10-2 

6.636 X 10-4 

3 

-1.508 X 10-1 

-1.594 X 10-2 

-4.597 X 10-4 

4 

9.830 X 10-2 

7.199 X 10-3 

-6.442 X 10-3 

5 

-6.946 X 10"2 

-4.046 X 10-3 

-4.103 X 10-5 

6 

5.047 X 10-2 

2.355 X 10-3 

-1.777 X 10-3 

7 

-3.761 X 10-2 

-1.515 X 10-3 

-9.176 X 10'5 

8 

2.814 X 10-2 

9.679 X 10-4 

-8.043 X 10-3 

9 

-2.118 X 10-2 

-6.537 X 10-4 


10 

1.582 X 10-2 

4.284 X 10-4 


11 

-1.174 X 10-2 

-2.881 X 10-4 


12 

8.571x10-3 

1.824 x 10-4 


13 

-6.148x10-3 

-1.142x10-4 


14 

4.286 X 10-2 

6.103x10-3 


15 

-2.893 X 10-3 

-2.628 X 10-3 

j 

16 

1.857 X 10-3 



17 

-1.117 X 10-3 



18 

6.060 X 10-4 



19 

-2.766 X 10-4 



20 

8.700 X 10-3 
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TABLE I.- RITZ COEFFICIENTS Cn FOR SEVERAL TYPICAL 
MINIMUM-MASS TRANSITION SECTIONS - Concluded 


(c) End conditions: Yq = 0.5; Xj = 1.5 




Cj^ for - 


n 

0) = -0.4 

0) = 0 

CD 

O 

II 

3 


m = 1.718 

m = 2.294 

m = 2.605 

1 

-5.936 X 10-1 

-8,310 X 10-2 

-2.800 X 10"2 

2 

1.696 X 10-1 

-1.850 X 10-3 

-2.491 X 10-3 

3 

-6.654 X 10-2 

-1.350 X 10-3 

-5.757 X 10-4 

4 

9.176 X 10-3 

-3.190 X 10-4 

-1.965 X 10"^ 

5 

1.271 X 10-2 

-1.661 X 10-4 

-8.378 X 10"® 

6 

-1.959 X 10“2 

-7.702 X 10“3 

-4.116 X 10"® 

7 

1.595 X 10-2 

-4.362 X 10"3 

-2.222 X 10"® 

8 

-9.262 X 10-3 

-2.571 X 10-3 

-1.272 X 10"® 

9 

2.277 X 10-3 

-1.595 X 10"3 

-7,473 X 10"® 

10 

2.321 X 10-3 

-1.098 X 10“3 


11 

-4.411 X10"2 

-6.224 X 10"® 


12 

4.183 X 10-3 

-6.929 X 10-® 


13 

-2.795 X 10"3 



14 

1.065 X 10“^ 



15 

2.350 XlO-*! 



16 

-9.189 X 10'^ 



17 

9.861 X 10-^ 



18 

-7.314 X 10"'^ 



19 

3.739 X 10-4 



20 

-1.227 X 10"4 
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TABLE II.- RITZ COEFFICIENTS Cn FOR MINIMUM-MASS 
TRANSITION SECTION WITH PRESSURE p = 0, 

Yq = 0.5, AND Xj = 1.5 


mac 

2-ny^2pFo 


2.309 


n 


1 

-7.731 X 10-2 

2 

-1.282 X 10-2 

3 

-4.310 X 10-2 

4 

-1.943 X 10"2 

5 

-1.029 X 10-2 

6 

-6.025 X 10-4 

7 

-3.777 X 10-4 

8 

-2.483 X 10-4 

9 

-1.688 X 10-4 

10 

-1.174 X 10-4 

11 

-8.283 X 10"^ 

12 

-5.870 X 10-5 

13 

-4.134 X 10-5 

14 

-2.844 X 10-5 

15 

-1.850 X 10-5 

16 

-1.030 X 10-5 
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TABLE m.- RITZ COEFFICIENTS Cn FOR MINIMUM-MASS 
TRANSITION SECTION WITH w = -1.0 AND = 1.5 

= 2.328 ' 

Yo = -2.784 


n 


1 

6.433 X 10-1 

2 

1.367 X 10-1 

3 

5.313 X 10-2 

4 

2.641 X 10"^ 

5 

1.499 X 10-2 

6 

9.216 X 10-2 

7 

5.970 X 10-2 

8 

3.998 X 10-2 

9 

2.730 X 10-2 

10 

1.879 X 10-2 

11 

1.288 X 10-2 

12 

8.654 X 10-^ 

13 

5.561 X lO-'l 

14 

3.255 X 10-‘l 

15 

1.931 X 10-4 

16 

1.103 X 10-4 

17 

5.488 X 10-2 1 
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TABLE IV.- RITZ COEFFICIENTS Cn FOR MINIMUM-MASS TRANSITION 
SECTIONS WITH STRESS RESULTANT CONSTRAINT 
= -0.4; zero-slope conditions; Xj = l.(^ 


(a) Design for Ng < 2.5; 
in = 1.565 


n 

Cn 

1 

-1.435 X 10“^ 

2 

-6.200 X 10-3 

3 

1.974 X 10"3 

4 

1.422 X 10-3 

5 

4.565 X 10-5 

6 

-3.718 X lO-*^ 

7 

-1.734 X lO-*^ 

8 

5.400 X 10~5 

9 

9.245 X 10-5 

10 

1.282 X 10-5 

11 

-3.666 X 10-5 

12 

-2.532 X 10-5 

13 

6.120 X 10-'^ 

14 

1.075 x 10-5 

15 

5.876 X 10-5 


(b) Design for no negative stresses; 
m = 1.277 


n 


1 

-3.478 X 10-1 

2 

1.030 X 10-1 

3 

-3.058 X 10-2 

4 

-7.938 X 10-^ 

5 

1.350 X 10-2 

6 

-1.650 X 10-2 

7 

1.443 X 10-2 

8 

-1.039 X 10-2 

9 

6.251 X 10-3 

10 

-3.011 X 10-3 

11 

9.632 X 10-4 

12 

2.698 X 10-5 

13 

-3.068 X 10-4 

14 

2.351 X 10-4 

15 

-8.651 X 10-5 





TABLE V.- NONDIMENSIONAL MASS m AND RITZ COEFFICIENTS Cn 
FOR MINIMUM-MASS PRESSURE -VESSEL HEADS WITH SLOPE AT 
HEAD-VESSEL JUNCTION EQUAL TO -0.5, 0, AND 0.5 



Cn for - 

n 

Slope = -0.5 

Slope = 0 

Slope = 0.5 


in = 0.7412 

m = 0.8503 

m = 1.249 


b = 0.5433 

b = 0.7739 

b = 0.9749 

0 

1.203 

1.387 

1.564 

1 

-1.828 X 10-1 

-3.350 X 10-1 

-2.334 X 10-1 

2 

-1.783 X 10-2 

2.038 X 10"^ 

-8.173 X 10-1 

3 

1.429 X 10-2 

-1.448 X 10-1 

7.024 X lO 'l 

4 

-3.483 X 10-2 

1.311 X 10-1 

1.111 

5 

-5.433 X 10-^ 

-5.914 X 10-2 

-2.414 

6 



1.520 X 10-1 

7 



1.894 

8 



-8.596 X 10"1 

9 



-9.918 X 10-2 
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